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Optimizing point-source parameters for scanning satellite 
surveys 
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ABSTRACT 

We describe a method for deriving the position and flux of point and compact sources 
observed by a scanning survey mission. Results from data simulated to test our method 
are presented, which demonstrate that at least a 10-fold improvement is achievable 
over that of extracting the image parameters, position and flux, from the equivalent 
data in the form of pixel maps. Our method achieves this improvement by analysing 
the original scan data and performing a combined, iterative solution for the image 
parameters. This approach allows for a full and detailed account of the point-spread 
function, or beam profile, of the instrument. Additionally, the positional information 
from different frequency channels may be combined to provide the flux-detection ac- 
curacy at each frequency for the same sky position. Ultimately, a final check and 
correction of the geometric calibration of the instrument may also be included. The 
Planck mission was used as the basis for our simulations, but our method will be 
beneficial for most scanning satellite missions, especially those with non-circularly 
symmetric point-spread functions. 

Key words: methods: data analysis - techniques: high angular resolution - tech- 
niques: photometric - surveys 



1 INTRODUCTION 

Survey missions designed primarily for observing the con- 
tinuum r adiation at one or more freq uencies (for e xample 
WMAP (|Bennett et al.lll997l l. Planck (|Tauberj|200"i )'). gen- 
erally produce as a by-product a point - or compact-source 
catal ogue (e.g. ILopez-Caniego et al.l 120071 ; IVielva et al.l 
l200&t ). containing positional and flux information on the 
detected sources. Generally these data have been, or are 
planned to be, extracted from the pixel maps on which the 
survey data are projected and collected. The process involves 
a point-source detection on those maps, followed by an es- 
timate of the position and flux from the pixels. There are 
two intrinsic problems here. Firstly, each pixel will represent 
a different collection of scan directions. Secondly, the digi- 
tisation on the pixel map will lead to signal distortion, as 
centres and sizes of map-pixels will never coincide accurately 
with the original size and shape of the measured samples. 
When the point-spread function (PSF) of the instrument is 
circular symmetric, these problems may be partly overcome. 
However, any asymmetry in the PSF will lead to a distortion 
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of the accumulated image that is difficult, if not impossible, 
to fully incorporate when deriving from the pixel maps the 
all-important image parameters: the position and flux of the 
source with their formal errors. It was also considered that 
pixelisation of the data could have a significant effect on the 
derived source parameters and their standard errors. 

Providing accurate and fully internally consistent po- 
sitional and flux information for point sources detected in 
a survey mission is crucially important for any subsequent 
use of those data. Cross identification followed by flux com- 
parisons for different wavelengths is only the most obvious 
application affected. Studies of the different spectral shapes 
among the detected sources will rely on the compatibility 
and accuracy of data obtained in different wavelengths to 
allow for reliable classification of (newly discovered) objects. 

We started the current study with the following hypoth- 
esis: The transit data used to build the point-source image 
contain together all information (image parameters) on the 
source, convolved for each transit with a PSF that is known 
or possible to reconstruct. Extracting the image parameters 
directly from the accumulated transit data should therefore 
make it possible to account directly for the applicable PSF 
for each transit. The question we then asked ourselves was 
the following: When extracting the image parameters directly 
from the transit data, will we see a significant improvement 
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in the accuracy and the statistics of the resulting image pa- 
rameters, compared to extraction from the accumulated data 
on a pixel map? To test this, one of us (ANM) prepared a 
software package with which the image parameter extrac- 
tion methods for pixel maps and for combined transit data 
could be tested and compared, while DLH prepared simu- 
lated data to be input t o the tests, following methods we 
used in an earlier study (|Harrison fc van Leeuwenll2005l ). 

The technique for extracting image parameters as pre- 
sented here, can be seen as a simplified version of the solu- 
tion that has been a pplied to the Hipparcos data reductions 
van Leeuwedl2007f ). Similar to those reductions, it can in 
principle also include a geometric calibration, which defines 
the relations between detector pointings at the sky and the 
satellite attitude. Similarly, as part of an iterative proce- 
dure, these methods can be used to prepare data as input 
for the calibration of the PSF. These kinds of dependencies 
are common for all-sky survey satellites, for which the data 
reduction is best treated as "self calibrating", to provide 
well defined internally rigid data systems. Such systems can 
then be calibrated as a whole to externally defined reference 
values. However, these calibration aspects are not treated 
in this paper, which concentrates on the extraction of the 
image parameters only. 

This paper has been organized as follows. Section [5] 
presents the mathematical background for the methods 
used, followed by a description of the data simulations in 
Section O The processing methods for the scan data to im- 
age parameters is presented in Section |4j and for the pixel 
maps in Section [5] The results of the application of these 
methods and a discussion of the results are presented in 
Section [6] and Section [T) 

We use ecliptic coordinates (A, 0) throughout the paper, 
because full-sky scanning strategies for satellites are gener- 
ally related to the ecliptic-plane position of the Sun, and 
tend to have symmetries with respect to the ecliptic plane. 



2 THE MATHEMATICAL MODEL 

2.1 The image response function 

The methods we describe here are applicable to sources for 
which the apparent diameter on the sky is small with re- 
spect to the PSF of the detector, but can in principle be 
extended to include sources for which the shape can, in first 
approximation, be represented by a known circular symmet- 
ric distribution. In both cases an image response function 
(IRF), R, can be created that predicts the image shape as 
the response of the instrument to an observation of a source. 

The IRF is described as an arbitrary function of the 
along- and across-scan coordinates, Aw = v t — v s and Tt- 
Here v s refers to the scan phase of the detector. The sub- 
script t refers to the coordinates of a target source, measured 
relative to the scan as described by the detector: 



R = R(Av,t) 



(1) 



The IRF incorporates the integration over fixed-length sam- 
ples along-scan, and is thus, for a point source, the actual in- 
strumental PSF convolved with a single-sample width block 
function. For the current application, the IRF is limited to 
a suitable radius, referred to as the image perimeter. The 



IRF is normalized in a relatively arbitrary way, for example 
by defining the peak to be at a level equal to 1.0. Similarly, 
the position of the peak can be defined as the centre of the 
IRF. Whatever way is chosen, it defines the observational 
response scale for the flux measurements, and the reference 
frame for the geometric calibration of the detector. The im- 
age response for an observed point source is simply given 
by the IRF times the intensity or flux of the source. A typ- 
ical transit for a source will consist of the responses for a 
sequence of equally spaced samples as obtained for a set of 
scan phases v s , roughly centred around the along-scan posi- 
tion of the source, Vt- This signal then includes the response 
from the source at the across-scan position of Tt- The scan 
circle itself is defined, for example, by the path described 
by the maximum possible response direction for a detector, 
and is usually close, but not equal, to a great circle. 



2.2 Projection effects 

The geometric calibration of an instrument such 
as the Planck satellite has been described in 
detail 



in 



Harrison fc van Leeuwen! (|2005l ) and 



lHarrison fc van Leeuwen! (|2006t ). where it is shown that an 



accuracy of about 0.01 times the FWHM (full-width half 
maximum) of the IRF can realistically be obtained. Thus, 
in the following discussion, uncertainties in the geometric 
calibration will in first instance be ignored, though we will 
come back to these very briefly in the final section of this 
paper. Similarly, the satellite attitude (in this case the 
position of the spin axis and the scan phase) can quite 
easily be determined to a relatively high accuracy, as this 
will be based on a star mapper and a high accuracy star 
catalogue. Thus, we will assume errors on the knowledge 
of the position of the detector as a function of time to 
be negligible with respect to other uncertainties, the most 
important of which will be the photon noise (or equivalent) 
on the original measurements. 

At each stage of the analysis we assume to have avail- 
able an estimated position and flux of the source we are 
investigating. In first instance this may simply come from 
the original identification of the source on a pixel map. The 
discrepancy between the actual and estimated positions and 
fluxes of the source will reflect in systematic differences be- 
tween the predicted and observed responses. These differ- 
ences together contain the information required for correct- 
ing the original to the best estimates of the source's position 
with respect to the scan (vt, Tt) and flux. 

The corrections to the estimated sky coordinates At , fit 
of the source are reflected in corrections to the Vt and Tt 
values for each scan, in a way that depends on the local 
orientation of the scan circle. The combination of scan circles 
with different local orientations allows for the reconstruction 
of the actual positional corrections. The local orientation of 
the scan circle is defined by the instantaneous position of 
the spin axis for the scan (X p ,fj p ), the opening-angle of the 
scan (the angle a between the detector and the spin axis) 
and the position of the source (\ t ,/3t)- The local inclination 
of the scan ip is now defined by: 



sin V = 



sin(At — A p ) cos /3 P 
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cos ip = 



sin 8 P — sin B t cos a 



(2) 



cos B t sin a 

with which the relations between Av, r and the corrections 
to the source-coordinates are expressed as: 

dvt — d\ t cos Bt cos ip + dBt sin ip, 

dr t = — d\ t cos ft sirup + dB t cosip. (3) 
2.3 The observations 

In the following we examine the data collected for a point 
source. We assume that we have an accurate predicted image 
shape, which is equivalent to the IRF for a flux 7": 



R P = / • R(v 



(4) 



This predicted image shape describes for each transit of the 
source the expected response as a function of the along- and 
across-scan positions of the source with respect to the scan 
circle. The signal is sampled at a series of discrete values j 
of v a , to provide observations: 



Oj — I ■ R(v t 



■T t ) +tj 



(5) 



where €j will generally be a function of the signal strength 
itself (such as Poisson noise). 

The predicted position, together with the attitude de- 
tails of the scans, defines for each scan predicted scan po- 
sition of the source, Vt and ft. Together with the predicted 
flux, this defines predicted observations for each scan, which 
can be compared with the observed counts: 

Oj = I ■ R(vt — v s j,f t ), (6) 



2.4 Correcting estimated values 

The differences between the predicted and observed instru- 
ment responses can be expressed as a first-order approxima- 
tion to corrections for vt, ft and I: 



+ 



dl ■ R(v t - v s j,f t ) 

- dR J ~ dR J 

1 ■ -Tz — dv t + 1 ■ tt— dn + ej , 

ov t art 



(7) 



where vt = vt + dvt, and similarly for the other variables. By 
combining Eq. and Eq. [7] the corrections to the assumed 
position and flux of the source, as derived from the observa- 
tions, are described as a set of linear equations, which can 
be solved by least squares: 



dl ■ R(v t 
'dR 



OR 



cos Bt cos ib — — — cos B t sin ib 

OVt OTt 



dXt 



dR . , OR 

- — sin ip + —— cos ip 

av t OT t 



dB t + €j 



(8) 



Potential weakness or ambiguity in the solution is resolved 
by the different orientations ip of the scan circles used, and 
the different values of the across-scan positions of the source 
as applicable to the different scans. Some iteration is still re- 
quired, as corrections to the three estimated image param- 
eters also affect the coefficients in Eq. [8] 



2.5 Extension to multiple frequencies 

The extension from a single to multiple frequency maps is 
simple, as long as the detectors for the different frequen- 
cies have been accurately calibrated with respect to each 
other. In that case, Eq. [8] is extended with a correction dl 
for each frequency (or possibly for each detector if there are 
suspected to be response differences between different detec- 
tors at the same frequency). The IRF has to be adapted for 
each detector, but the corrections dAt and dBt are the same 
for all frequencies. Extension to multiple frequencies could 
bring improvement in positional accuracy, and so indirectly 
also to flux estimates. 



2.6 Solving for a background gradient 

A linear 2-dimensional gradient in the background can be 
resolved through a minor extension of Eq. [5] The reason 
for solving for a possible linear background gradient is its 
effect on the estimated position of the source, which will get 
biased by an unsolved background gradient. The background 
gradient is expressed in offsets for the sky-coordinates of the 
source: 



013 dB 
o{\cosB) op 



(9) 



The quantities dA cos 8 and dB are directly related to the 
position of the detector relative to the source through the 
orientation angle ip of the scan as defined in Eq. 



dAcos/3 = cos ip(v S:J — v t ) — sin ip(jt) 
dB = sin ip(v s j — v t ) + cos ip(r t ) 



(10) 



Substituting Eq.[TD]into Eq.[5]describes the background con- 
tribution as a function of the source observation coordinates. 
As such it can be added to Eq.[8]as an additional set of un- 
knowns to be solved. Some care is required here though, 
as adding more degrees of freedom has a potential for de- 
stabilizing the solution, in particular when relatively few 
independent observations are available. Independent here 
refers primarily to the distribution of scan direction. The 
standard errors on the background coefficients as derived 
from the least squares solution, as well as the standard de- 
viation of that solution, can be used as criteria to decide 
whether adding a background gradient as additional param- 
eters to the solution is at all beneficial. 



3 THE DATA SIMULATIONS 
3.1 Scanning law 

Full-sky scanning missions with severe constrai nts on instru- 
ment temperature co ntrol, such as Hipparco s |van Leeuwenl 
Planck, Gaia (|Lindegren et al.ll2008l ). have scanning 
laws implemented that maintain a fixed angle between the 
spin axis of the satellite and the direction to the Sun. This 
angle is referred to as the solar aspect angle £. The full-sky 
scanning is obtained by two motions: a precession of the 
spin axis around the (nominal) direction of the Sun, and a 
rotational motion of the instrument around the spin axis. 
This type of scanning strategy has been referred to as "cy- 
cloidal" . The limitation of such a scan is the coverage close 
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to the ecliptic plane, for which the spread in inclination an- 
gles for the scan coverage is limited by £. If the inclination 
of a scan circle with respect to the ecliptic plane is given 
by cp (with if = ±7r/2 for a scan circle crossing the ecliptic 
plane at right angles) then 



| cos <p\ ^ 



sin £ 



(11) 



where a is the opening-angle of the detector as defined above 
in Section r2.2l Thus, for a small solar aspect angle, the incli- 
nations of the scan circles remain all close to perpendicular 
to the ecliptic plane. This is a relatively unfavourable sit- 
uation for point- and compact-source analysis as the easier 
to extract positional information is in the along-scan direc- 
tion. This is therefore a good configuration to see how the 
current analysis will perform under difficult conditions. For 
missions such as Hipparcos and Gaia, which aim specifically 
at the astrometry of point sources, the Solar-aspect angle 
has be en maximized, at values of 43 and 45 degr ees respec- 
tively (|van Leeuwenlfl997l ; iLindegren et alj|200ah . while for 
Planck i t has been chosen at a r elatively small value of 7.5 
degrees (|Dupac fc Tauberl |2005j ). The Planck mission has 
therefore been chosen as an example case for our experi- 
ments, as it provides a wide range of conditions depending 
on the ecliptic latitude of the source. The precession period 
for the spin axis of the Planck satellite is scheduled to be 6 

months. 

The Planck scanning (|Dupac fc Tauberl I2005T) is de- 
scribed not as a continuous function, as for Hipparcos and 
Gaia, but step-wise. This means that the spin axis remains 
at a "fixed" position over a period of, for example, 1 hour, 
during which the satellite describes 60 full circles. The nom- 
inal scan velocity is 6 degrees per second, and every part of 
the scan is examined by 60 successive scans. The accumu- 
lated data from the scan circles are referred to as a "ring" . 
Due to the time required for repositioning the spin axis, 
there are effectively only about 55 useable scans per point- 
ing. Scans in successive pointings will have a maximum sep- 
aration of about 2.5 ± 0.65 arcmin. The mean value is the 
mean motion of the Sun on the ecliptic, the range is the re- 
sult of precession of the spin axis around the Sun, as a result 
of which the scan density will be higher when the spin axis 
moves in opposite direction of the Sun and vice versa. For 
our experiments we have used one year of "mission data", 
providing a generally homogeneous scan coverage of the sky. 



3.2 The detectors 

For testing the methods described in Section [2] we used a 
"worst-case" scenario by selecting the 30 GHz detectors, 
which have the largest beam size (FWHM= 33 arcmin) and 
the largest predicted beam ellipticity (1.4). For this detector 
we adapted an opening-angle of 85 degrees (this is for test- 
ing purposes only and doesn't represent the actual detector 
position in the Planck focal plane). At the sampling rate 
of 32.5 Hz, the samples are spaced 5.4 arcmin apart, and a 
point-source image is effectively covered by 6 to 8 samples. 

The IRF (see section [2| has been represented by a 2- 
dimensional Gaussian distribution: 



R(v, t) = exp 



Av 2 
' 2a?, 



T 

2a* 



(12) 



where a = 0.5/^2ln2 w 0.425 FWHM. The shape of the 
IRF has only a very minor effect on the results presented 
in this study. The main effect on the positional accura- 
cies comes from the image "width" and the steepness of its 
slopes, and a Gaussian profile forms a good first approxima- 
tion for a typical beam profile. 



3.3 Samples and sample noise 

Two distributions of samples have been considered, one fol- 
lowing a realistic scan-velocity, which will generally cre- 
ate an uneven distribution of samples over an image, and 
a "smooth" sampling, with all samples spread evenly over 
the image. The latter was finally preferred for testing, as it 
should provide an simple-to-interpret set of results, not un- 
necessarily complicating the possible detection of systematic 
errors that may possibly be inherent to the methods tested. 
Thus, at each pointing direction a set of evenly spaced sam- 
ples was created that represented the 60 scans made at that 
pointing. Gaussian noise was added to these samples at a 
level according to a specified signal to noise ratio, which 
applied to the peak of the image. The contributing scans 
for each source were determined by the adopted scanning 
law, which defines the orientation and position of each scan 
relative to that of the source. For each set of contributing 
transits of a given source, 100 different noise realizations 
were made. 



4 THE PROCESSING OF THE SCAN DATA 

The available simulated data for each source were inves- 
tigated following the methods described in Section 12.41 
through weighted least-squares solutions, providing esti- 
mated values for the position and flux (the image param- 
eters), and the standard deviations on those values. These 
parameter values were compared with the input data, and 
the differences compared with the estimated standard er- 
rors. The weights used in the solutions were based on the 
standard deviation of the Gaussian noise used in the data 
simulations. The mean values of the image parameters and 
their standard errors were determined as based on the 100 
independent noise realizations. This allows in principle for 
bias detection at 10 per cent of the standard deviation. For 
a S/N = 10, typical values of just under 2 arcsec were found 
for the positional standard error, with a sigma for the beam 
of 841 arcsec. For the flux, the standard error was approx- 
imately 0.15 per cent. The observed dispersion in the mea- 
sured values was in good agreement with the derived stan- 
dard errors for the solutions. More details on results are 
provided in Section [6] 



5 IMAGE PARAMETERS FROM PIXEL MAPS 

As part of the experiment, the samples as obtained in the 
scans were projected on a map of square pixels. To avoid 
complications caused by distorted pixels, each source was 
projected as if it was situated at the ecliptic plane, where 
the map approaches a flat patch. The pixel size was chosen 
at 6 arcmin, which is compatible with the sample size on the 
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scans (5.4 arcmin) and the maximum distance between suc- 
cessive scans (3.15 arcmin), and guarantees that at least two 
samples are contributing to each pixel. The 6 arcmin pixels 
correspond to what is r eferred to in the H EALPbfl software 
as NSIDE = 2 9 = 512 (|G6rski et alj|2005h . At the next res- 
olution level, the pixel size would be 3 arcmin, leaving some 
pixels empty. 

Two options were explored, one in which the source 
took a random position within a pixel, and a second where 
the source was central in a pixel. This allowed for some dis- 
tinction of how pixelization may affect the extracted image 
parameters. 

After defining the position of the source with respect 
to the map pixels, the samples from the accumulated scans 
were projected on the relevant pixels, where each sample 
was assigned in full to the nearest pixel, as is standard 
procedure for creatin g CMB pixel maps (see for example 
lAshdown et al.ll2007f ). At each pixel, the response was sim- 
ply taken as the mean of all contributing samples. 

The image-parameter extraction for the pixel maps fol- 
lows the same principles as for the scan data. We assumed in 
all cases a circular symmetric 2-dimensional Gaussian pro- 
file for the image, with the width defined by the mean a of 
the actual beam. In case of an asymmetric beam profile this 
will be somewhat inaccurate, but in that situation the image 
profile on the map is always poorly defined due to the differ- 
ent contributions to individual pixels. Thus, we start with 
a 2-dimensional Gaussian profile for the IRF as in Eq. 1121 
but now with respect to pixel coordinates (x, y) and a first 
assumption of the source position (xo,yo)' 



R(x,y) = exp 



(x-x f (y-y ) 2 



2a 2 



2a 2 



(13) 



The predicted response for a source in pixel (x,y) is simply 
given by: 



5 x ,y = / ■ exp 



(x-x ) 2 (y-yo) 2 



2a 2 



2a 2 



(14) 



Corrections to the estimated image parameters (xo,yo,I) 
are obtained in the same way as was done for the scan data, 
by comparing the predicted with the observed responses: 

OR - OR 

Ox, y - o x , y = R(x, y) dl + I —-dx + I -7—dy. (15) 

ox ay 

These equations are solved by least squares. Iterations are 
required as corrections for the position estimate affect the 
predicted responses in a non-linear way. The final iteration, 
based on a convergence criterion, provides a set of instru- 
ment parameters and their standard errors. 



6 RESULTS AND COMPARISONS 

In the comparisons of results, the following performance cri- 
teria were looked at: 

• Accuracy of the reproduction of the source position; 

• Accuracy of the predicted standard error on the source 
position; 
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• Accuracy of the reproduction of the flux; 

• Accuracy of the predicted standard error on the flux. 

Together these criteria define both the reliability and ro- 
bustness of the methods that have been applied. 

6.1 Distribution of test sources 

Tests were carried out with sources distributed evenly over 
all ecliptic longitudes at 1 degree separation and at a few se- 
lected fixed ecliptic latitudes. The scan coverage of a source 
is a strong function of both ecliptic latitude and longitude, 
and any algorithm applied to source-parameter extraction 
will need to be able to provide reliable information under 
all these conditions. This is particularly the case for a rel- 
atively short mission such as Planck, where the mission is 
expected to last less than 2 years. The local variations in 
the distribution of scan directions also become more pro- 
nounced for a Planck-like scan for which the solar aspect 
angle is small. 

6.2 The test data sets 

The test data sets we created aimed at testing the methods 
under increasingly more complex and difficult conditions. 
The first data set, described in Section \6. 31 used a circular- 
symmetric beam and S/N = 10. Map pixelization was done 
by randomly positioning of the source with respect to the 
pixels on the "map" . Data were obtained at ecliptic latitudes 
of and 60 degrees. 

In Section [6.41 the pixel positions were forced such that 
the source would be positioned exactly at the centre of a 
pixel. Here data were only obtained a latitude of degrees. 

Finally in Section 16.51 the most extensive tests are de- 
scribed. These covered 4 latitude (0, 30, 60 and 89 degrees), 
4 ellipticities (1.0, 1.2, 1.4 and 1.6) and 7 signal-to-noise ra- 
tios (S/N = 0.1, 0.2, 0.5, 1.0, 3.0, 5.0, 10.0), equivalent to an 
intensity range from 0.01 to 100. There was only one lon- 
gitude used, A = 240 degrees, which had shown in circular- 
beam tests to be optimal for the map-making results, inde- 
pendent of pixel positions. 

6.3 Circular-symmetric beam 

This data set provides what should be the easiest conditions 
for the pixel-map extractions, and can show the minimum 
difference in performance between the two methods. 

6.3.1 The scan-based analysis 

All results for the scan-data analysis are internally consis- 
tent and at, or very close to, the statistical limit. The mean 
standard errors over the 100 noise realizations as a function 
of ecliptic longitude display the variations in density of scan 
coverage (Fig. [T]). One of these variations is caused by the 
smaller than average displacements between scans when the 
spin-axis movements are opposite to the apparent movement 
of the Sun (causing the denser scan at around 55 and 235 
degrees) . As was stated already in Section 13.11 when those 
movements are in the same direction, the scan is less dense 
and the standard errors are larger (at longitudes of 145 and 
325 degrees). Also visible are a couple of data points affected 
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A (degr) 

Figure 1. Standard errors for scan-based solutions as a function 
of ecliptic longitude, in both graphs for latitudes of (top) and 
60 (bottom) degrees. The lower graph shows the results in longi- 
tude (crosses) and latitude (dots), which are observed to be very 
similar. This could be expected for a circular symmetric beam. 
The few points showing in each graph higher than average accu- 
racies are due to an overlap at the closure of the full scan. The 
top graph shows the same for the flux determinations, where the 
errors need to be compared with a peak flux of 10.0. 

by the "full closure" of the scan, causing effectively a few 
additional scans covering longitudes at 20 and 200 degrees. 
Very similar curves are observed for the standard errors on 
the flux, as is shown in the top graph of Fig. [T] A a value 
of 0.02 is equivalent to a 0.2 per cent error on the flux mea- 
surement. The next test is to see how the mean standard 
errors for each longitude compare with the observed stan- 
dard deviations, as based on the analysis of the 100 noise 
realizations (Fig. [2}. The spread of the data points is close 
to what could be expected for standard deviations based on 
100 independent measurements. We observe a small increase 
of the standard deviations for the lowest fluxes. The same 
is observed for the errors on the flux measurements. There 
the increase accompanies the development of a positive bias 
in the flux estimates. 

The final comparison is between the input and observed 
values for the positions and fluxes. Values have been ob- 
tained while averaging over all longitudes. Small and in- 
significant differences are seen depending on whether the ob- 
served standard deviations or the mean standard errors are 
used as weights. Given the close agreement between these 
quantities, this is no surprise. The errors on the mean val- 
ues per longitude have been taken as 0.1 times the mean 
standard error, considering that 100 observations were used 




1 2 



a arcsec (Mean standard ernor) 

Figure 2. Comparison between the mean standard errors and 
observed standard deviations for 100 noise realisations for the 
scan-based solutions. The data in longitude (crosses) and latitude 
(dots) are shown for both the latitude and 60 simulations. 



Table 1. Reproduction of input data with scan-based solutions 



Coord. 


Resid. 
(3 = 


uwsd 


Rcsid. 
(3 = 60 


uwsd 


<5A cos (3 

5(3 
r-io 

0.01 


0.003 ±0.012 
0.012 ±0.012 
-0.003 ± 0.010 


1.024 
0.932 
0.943 


0.011 ± 0.008 
-0.010 ±0.008 
-0.009 ± 0.007 


1.009 
1.003 
0.972 



uwsd:unit-weight standard deviation 



to calculate each of the means. The results are summarized 
in Table [TJ which shows that there are no biases observed in 
the scan-based analysis results for either positions or fluxes. 

6.3.2 The pixel-map based analysis 

The results for the pixel-map based analysis are less accurate 
or predictable. We observe a poor relation between the stan- 
dard errors of the solutions and the actual observed standard 
deviations, as shown in Fig. [3] The comparison shows that 
the standard errors derived from the least-squares solutions 
generally provide in this situation a very poor estimate of 
the actual errors on the measurements. A comparison with 
Fig. [2] shows that for the analysis of the pixel-map data the 
positional errors are mostly larger, by up to a factor ten, 
than the scan-based analysis of the same data. The main 
problems with the analysis of the pixel-map data are caused 
by poor scan coverage. Each source on the ecliptic plane will 
generally receive two scans over a 12 months period. At the 
ecliptic plane these scans will have a difference in inclina- 
tion by between and 15 degrees (twice the solar aspect 
angle of the spin axis). This is referred to as the intercept 
angle. For values close to zero, the standard deviation for 
the solution tends to be large, while at maximum value the 
pixel-map based results become close to those of the scan- 
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Figure 4. Lower graphs: the observed standard deviations for pixel-map based solutions at latitude (left) and latitude 60 (right). 
Upper graphs: the scan intercept angle (see text) as a function of ecliptic longitude. Dots and crosses refer to measurements in latitude 
and longitude respectively. 



Table 2. Reproduction of input data with pixel-map based solu- 
tions 



Coord. 


Resid. 
/3 = 


uwsd 


Resid. 
P = 60 


uwsd 


<5A cos p 
5/3 

J-10 
0.01 


-0.019 ±0.021 
0.002 ± 0.025 
-7.4 ±0.18 


0.957 
1.412 
0.943 


-0.115 ± 0.011 
0.016 ± 0.012 
-7.4 ±0.093 


1.361 
2.393 
2.298 



uwsd:unit-weight standard deviation 



based analysis (Fig. [4j. The same effects are observed in 
the errors for the flux measurements. However, while there 
is no bias in the positional measurements observed, there 
is a significant flux loss (see Table [2}. A simple compari- 
son between the data in Tables [T] and [2] as well as between 
Figures [5] and [3] clearly demonstrates the performance ad- 
vantages of the scan-based analysis, providing reliable and 
generally smaller standard errors and a bias-free flux esti- 
mate, with no noticeable dependence on the intercept angle 
(strictly speaking, an intercept angle can only be defined 
when two sequences of related scans are present). 



6.4 Circular symmetric beam, source position 
fixed at pixel centre 

In order to investigate the origin of the rather large stan- 
dard deviations observed in the map-based solutions, a set 
of solutions at latitude degrees was made such that the 
image was always positioned at the centre of a pixel on the 
map. The standard deviations of the positional fits are now 
closer to what is found for the scan data (Fig. [SJ), though 
the noise level on the map data is on average higher and 
repeatability of those data is a lot poorer. The relation be- 
tween those standard deviations and the standard errors is 
in addition still very poor too. The mean standard error 
estimates following from the least squares solutions of the 
map data are systematically overestimated and generally 
very noisy (Fig. [6} . Thus, even in these, what may appear 
to be relatively favourable, but still unrealistic, conditions, 
the pixel-map derived positional information is still poor, in 
the prediction of the errors, while the errors themselves are 
poor estimates of the actual noise on the data. This can be 
observed in Fig. [7] where the offsets for the recovered po- 
sitions in longitude and latitude are shown as a function of 
longitude of the source. For comparison, the results for the 
scan-based analysis are also shown. An example of how the 
distribution of pixels affects the data accumulated from the 
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Figure 3. Comparison between the mean standard errors and 
observed standard deviations for the pixel-map based solutions. 
The diagonal line shows the one-to-one relation. Dots and crosses 
refer to measurements in latitude and longitude respectively. 
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Figure 5. Positional standard deviations over 100 noise realisa- 
tions, as a function of ecliptic longitude, for sources at ecliptic 
latitude zero. The broad distribution of dots are the results for 
pixel-map based positions, where the source was always forced 
to be at the centre of a pixel. The narrow distribution of crosses 
shows the same data analysed directly from the scans. 
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Figure 6. Observed positional standard deviations compared 
with mean standard errors from the image-parameter solutions. 



Table 3. Reproduction of input data with pixel-map based solu- 
tions for fixed source positions 



Coord. 


uwsd 
Pixel- 


sd 
map 


uwsd sd 
Scans 


Units 


<5A cos /3 
6(3 

1-10 
0.01 


1.106 
2.518 
1.547 


4.92 
10.43 
5.50 


1.024 0.237 
0.932 0.213 
0.972 0.185 


arcsec 
arcsec 



scans is shown in Fig. [8] A summary of these results as well 
as a comparison with the scan-based analysis is presented 
in Table Here both the actual standard deviation and the 
unit- weight standard deviation are shown for the same data 
analysed from pixel maps and directly from scan data. The 
conclusion can be drawn that the scan-based analysis is per- 
forming more than an order of magnitude better than the 
pixel-map based analysis, in the reconstruction of all image 
parameters as well as in the provision of standard errors. 

Also in flux reconstruction the results for the pixel-map 
derived solutions remain poor, as is shown in Fig. [7] There 
is still a systematic flux loss of on average about 7 per cent, 
and a noise that is in particular high for small intercept 
angles. None of these features shows in the analysis directly 
made from the scan data, the results of which are also shown 
in Fig. [7] Conclusions we can draw from these experiments 
are that the position of the source centre with respect to the 
pixels on the pixel map adds a significant source of noise, 
in particular for small intercept angles. But even when this 
"positional freedom" is removed, and the sources are artifi- 
cially placed at pixel centres, the results remain significantly 
inferior to those obtained directly from the scan data, in re- 
production of the image parameters (with much higher noise 
and a flux loss for the map-based analysis) and in accuracy 
of the standard error estimates. Considering that we are still 
dealing here with a circular-symmetric beam, it appears that 
the actual pixelisation process is a major source of additional 
noise in the source-parameter determination. 
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Figure 9. Standard deviations of positions for different ellipticity, based on the direct analysis of the scan data, as a function of source 
intensity. Data in longitude direction are shown as crosses, in latitude direction as dots. The differences between the longitude and 
latitude results for higher ellipticity values reflect the average beam shape. In this example the source was positioned at 30 degrees 
latitude. 



6.5 Varying ellipticity, flux and latitude 

For the final test runs, on the effects of ellipticity on the 
image-parameter reconstruction, an ecliptic longitude of 
240 degrees on the sky was chosen for which the pixel-map 
based analysis was giving the relatively best results. Sources 
were placed at latitudes of zero, 30, 60 and 89 degrees. The 
fluxes of sources at those positions were varied over a range 
of S/N values from 0.1 to 10, equivalent to a range of source 
flux from 0.01 to 100. The ellipticity, defined as the ratio of 
the FWHM along scan over the FWHM across scan, was 
tested for values of 1.0, 1.2, 1.4 and 1.6, where a value 
of 1.4 is typical for the lowest freq uency detectors on the 
Planck satellite (|Sandri et alj [20051 ). Each point in this 3- 
dimensional data grid was simulated with 100 independent 
noise realizations, The simulated data were analysed from 
the scans and projected onto pixel maps. 

Figure [9] shows an example of the results from the scan- 
based analysis, in this case for a source at latitude 30 de- 
grees. Results from the scan-based analysis for other lat- 
itudes were very similar and mainly reflected difference in 
total scan coverage. As can be observed from Fig.|5J the stan- 
dard deviations of the fitted positions as obtained from the 
different noise realizations behaves fully in the expected way, 



i.e. they are proportional to the square root of the flux. The 
only additional feature visible is a split between the stan- 
dard deviations for the longitude and latitude directions for 
the high-ellipticity beam profiles, which simply reflects the 
preferential direction of scans being close to perpendicular 
to the ecliptic plane. The effective beam width in longitude 
is therefore narrower than in latitude. 

Equally important for the analysis is the reliability of 
the standard errors. A comparison between the standard er- 
rors and the observed standard deviations for all test cases 
as obtained from the scan-based analysis is shown in Fig. 1101 
No dependencies on latitude or ellipticity are observed, and 
there is only a small increase of the observed standard devi- 
ations compared to the mean standard errors for the faintest 
sources. 

The errors on the reconstructed flux values are shown 
in Fig. 1111 A small bias is observed for the lowest two flux 
levels, which doesn't show a dependence on either latitude 
of the source or ellipticity of the beam. No significant bias 
is observed for sources at flux level 1 and above, i.e. those 
with S/NJS 1 on the original scan data. The errors on the 
fluxes are equivalent to an accuracy of order 0.01 magn. or 
better for all sources with a S/N^ 1 in the scan data. 

Exactly the same data has been analysed after projec- 
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Figure 12. The same data as shown in Fig. [9] but now analysed from the pixel map. 



tion on a pixel map as described earlier in this paper. No 
special provisions were made about placement of the source 
position with respect to the pixels, but the longitude of the 
sources had been chosen such that what appeared to be the 
most favourable conditions for the pixel-map based solutions 
would be obtained. At almost any other longitude the tests 
described in Section 16.31 showed results that are (consider- 
ably) worse. In contrast, a similar variation in longitude was 
shown to have only very small effects (only due to total cov- 
erage variations) on the direct analysis of the scan data. 

Figure [12] shows the equivalent of Fig. |9j but now for 
the pixel-map based analysis. A couple of differences can 
be noticed. For the brightest images the sensitivity to the 
ellipticity of the beam has an increasingly damaging effect, 
which is what could be expected, as the image-parameter 
fits for brighter images become increasingly more sensitive 
to an accurate beam representation. For all sources the er- 
rors are somewhat larger than for the scan-based analysis, 
but this is in particular the case for the fainter ones. It should 
be noted too that in the pixel-map based analysis between 
4 and 8 per cent of the parameter fits on the faintest im- 
ages were unsuccessful (i.e. didn't converge). No such fail- 
ures were found in the scan-based analysis. 

As was observed for the circular beam analysis, the 
standard errors obtained from the pixel-map based solutions 
are inaccurate, and tend to be, at this longitude, an overesti- 



mate. This is seen here in Fig. 1131 which, however, also shows 
that for the faintest images the actual standard deviations 
can be larger than the mean standard errors obtained from 
the solutions, similar to what is observed for the scan-based 
analysis. 

Figure [14] shows the reconstruction of the flux values 
for the pixel-map based analysis, and should be compared 
with Fig. [TT]of the scan-based analysis. As was observed in 
Section [6. 31 the reconstructed fluxes in the pixel-map based 
analysis are biased, systematically providing too low esti- 
mates of the actual flux. This becomes worse in the case 
of an elliptical beam, and in particular close to the ecliptic 
plane, where the image on the map is constructed from only 
scans at a narrow range of intercept angles. In magnitude 
terms, the bias would be of order +0.1 for an ellipticity of 
1.4 to 1.6. 

In conclusion, the pixel-map based analysis for even the 
most favourable map data is still significantly less accurate 
than the scan-based analysis for any of the data. This applies 
to all aspects of the image-parameter extraction: positions, 
fluxes and the standard errors on the estimated parameters. 



7 DISCUSSION AND CONCLUSIONS 

Next to the frequency maps, the compact source catalogue 
provides a lasting heritage of a survey mission such as 
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Figure 8. Examples of a pixel grid coverage of scan observations, 
showing the coverage variations due to accidental positioning of 
the grid with respect to the source position. The central dot rep- 
resents the source position. The connected sequences of small dots 
represent samples on successive scans. 
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Figure 10. A comparison between measured standard deviations 
and mean standard errors for the scan-based solutions for all test 
cases. Only for the faintest sources (at a S/N on the scan of 
0.01) is there a small increase in the observed with respect to the 
predicted noise level. The crosses and dots refer to the data for 
longitude and latitude respectively. 
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Figure 7. Lower three graphs, from top to bottom: flux measure- 
ments, residuals in latitude and residuals in longitude. The wide 
spread of points in each graph originate from the pixel-map based 
analysis, the narrow distributions from the scan-based analysis. 
Upper graph: the intercept angle of the scans contributing to the 
data. All data are shown as a function of ecliptic longitude and 
concerns sources at ecliptic latitude zero. 



Planck. It should therefore be prepared with the utmost care 
and to the best capacity of the actual data, so that it can be 
reliably compared with survey data and source catalogues 
in other frequency domains. 

It appears that the main source of noise in the pixel- 
map based extraction of image parameters originates from 
the manner in which the data are projected on the pixel 
map. However, any more sophisticated method, considering 



Figure 11. Recovery of the input flux in the scan-based solution 
as a function of the source flux on a single scan. All 16 test cases 
(4 latitudes times 4 ellipticity values) are shown for each flux. 
Multiplication of the vertical scale by 2.5 gives the equivalent 
errors in magnitudes. 
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Figure 13. A comparison between mean standard errors and 
observed standard deviations for the pixel-map based analysis. 
Data are shown for all test cases. The equivalent data for the 
scan-based analysis is shown in Fig. 1101 The crosses and dots 
refer to the data for longitude and latitude respectively. 



Figure 14. The equivalent of Fig. 1111 but now for the pixel- 
map based analysis. The spread in data points for each flux value 
is closely related to the ellipticity (indicated by colour, from low, 
red, to high, blue), with a secondary dependence on latitude, such 
that the worst bias is, as could be expected, found for the highest 
ellipticity, while sources at higher latitude are less affected than 
those at lower latitudes. 



for example a spread of data over several pixels, will be com- 
plicated and confuse the statistical properties of the data. 

We have shown that there is a significant advantage in 
extracting source parameters from a full-sky scanning survey 
mission directly from the scans rather than from the data as 
projected on a pixel map. The methods themselves are sim- 
ple to implement, are stable, and provide accurate results for 
positional and flux information on the sources as well as re- 
liable standard errors on these parameters. The methods as 
described here are easily extendible to multiple frequencies 
and can also take into account background gradients. These 
are trivial extensions of the implementation we have tested 
here. As a further extension, but somewhat more compli- 
cated and requiring further tests, an implementation could 
also include a final adjustment of what is generally referred 
to as the instrument parameters, describing the relative po- 



sitions and flux-responses of a set of different detectors in 
the focal plane. 

Extending the analysis of the image parameters across 
different frequency maps has the added advantage of pro- 
ducing source information that is always consistent across 
different maps, and could even take fully into account situ- 
ations where two sources are resolved on a high-resolution 
map, and appear merged on a low resolution map, providing 
well-defined separate fluxes on all maps. 

A further advantage of the scan-based over the map- 
based analysis is the ability of the former to determine the 
source flux over a range of epochs, determined by the fre- 
quency the source is transitted rather than the frequency at 
which the maps may be made. 

A good strategy for the creation of a complete and re- 
liable source catalogue for a survey mission could be the 
following: 

(i) After a proper geometric calibration of the instrument, 
prepare the frequency maps from the scan data; 

(ii) In the maps, at all frequencies, identify any feature 
that might be due to a compact source. This catalogue may 
be overcomplete, in other words may contain spurious de- 
tections, since these will be resolved and removed by the 
scan-based analysis; 

(iii) Extract a list of potential source positions for all 
maps together; 

(iv) Extract for each potential source the scan data at 
each frequency that would be within the beam-range of the 
source; 

(v) Apply the scan-based solution on the accumulated 
data for all frequency maps together, resolving any merged 
images where possible on the basis of higher resolution data; 

(vi) Optionally, examine residuals for individual detectors 
on systematics that can reflect small corrections to the ge- 
ometric and response calibrations; if these are found to be 
significant, apply and rebuild the pixel maps, then repeat 
points (ii) to (v). 

What will be obtained this way is a compact source cata- 
logue that is reliable in positional and flux data, as well as 
the errors assigned to those data, and which is in addition 
fully consistent over all frequencies and very likely to be as 
complete as it can possibly be within the detection limits. In 
other words, a compact source catalogue that can have last- 
ing value. It is also a catalogue that will automatically be 
free from any ambiguities between the results for different 
frequency maps. 



ACKNOWLEDGMENTS 

Some of th e results in this pap er have been derived using the 
HEALPix (|G6rski et al.ll2005T l package. ANM is grateful for 
the support provided by a Marshall Scholarship and an NSF 
Graduate Research Fellowship during the preparation of this 
work 



REFERENCES 

Ashdown M. A. J. 
A&A, 467, 761 



Baccigalupi C, Balbi A., et al. 2007, 



Point-source parameters for scanning surveys 



Bennett C. L., Halpern M., Hinshaw G., Jarosik N., Limon 
M., Mather J., Meyer S. S., Page L., Spergel D. N., Tucker 
G., Wilkinson D. T., Wollack E., Wright E. L., 1997, in 
Bulletin of the American Astronomical Society Vol. 29 
of Bulletin of the American Astronomical Society, The 
Microwave Anisotropy Probe (MAP) Mission, p. 1353 

Dupac X., Tauber J., 2005, A&A, 430, 363 

Gorski K. M., Hivon E., Banday A. J., Wandelt B. D., 
Hansen F. K., Reinecke M., Bartelmann M., 2005, ApJ, 
622, 759 

Harrison D. L., van Leeuwen F., 2005, MNRAS, 364, 1441 
Harrison D. L., van Leeuwen F., 2006, MNRAS, 365, 699 
Lindegren L., Babusiaux C, Bailer- Jones C, Bastian U., 
Brown A. G. A., Cropper M., H0g E., Jordi C, Katz 
D., van Leeuwen F., Luri X., Mignard F., de Bruijne 
J. H. J., Prusti T., 2008, in IAU Symposium Vol. 248 of 
IAU Symposium, The Gaia mission: science, organization 
and present status, pp 217-223 
Lopez-Caniego M., Gonzalez-Nuevo J., Herranz D., Mas- 
sardi M., Sanz J. L., De Zotti G., Toffolatti L., Argiieso 
F., 2007, ApJS, 170, 108 
Sandri M., Villa F., Burigana C, Mennella A., Nesti R., 
Mandolesi N., Bersanelli M., 2005, in Mazuray L., Wart- 
mann R., eds, Society of Photo-Optical Instrumentation 
Engineers (SPIE) Conference Series Vol. 5962 of Society of 
Photo-Optical Instrumentation Engineers (SPIE) Confer- 
ence Series, Straylight analysis and minimization strategy 
in PLANCK low frequency instrument, pp 575-586 
Tauber J. A., 2004, Advances in Space Research, 34, 491 
van Leeuwen F., 1997, Space Sci. Rev., 81, 201 
van Leeuwen F., 2007, Hipparcos, the new reduction of the 

raw data. Springer, Dordrecht 
Vielva P., Martinez-Gonzalez E., Gallegos J. E., Toffolatti 
L., Sanz J. L., 2003, MNRAS, 344, 89 



